home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / SSIDI.z / SSIDI
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. SSSSSSSSIIIIDDDDIIII((((3333FFFF))))                                                            SSSSSSSSIIIIDDDDIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SSIDI   - SSIDI computes the determinant, inertia and inverse of a real
  10.      symmetric matrix using the factors from SSIFA.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.       SUBROUTINE SSIDI(A,LDA,N,KPVT,DET,INERT,WORK,JOB)
  15.  
  16. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  17.      On Entry
  18.  
  19.      AAAA REAL(LDA,N)
  20.         the output from SSIFA.
  21.  
  22.      LLLLDDDDAAAA INTEGER
  23.         the leading dimension of the array A.
  24.  
  25.      NNNN INTEGER
  26.         the order of the matrix A.
  27.  
  28.      KKKKPPPPVVVVTTTT INTEGER(N)
  29.         the pivot vector from SSIFA.
  30.  
  31.      WWWWOOOORRRRKKKK REAL(N)
  32.         work vector.  Contents destroyed.
  33.  
  34.      JJJJOOOOBBBB INTEGER
  35.         JOB has the decimal expansion  ABC  where
  36.         If  C .NE. 0, the inverse is computed,
  37.         If  B .NE. 0, the determinant is computed,
  38.         If  A .NE. 0, the inertia is computed.
  39.         For example, JOB = 111  gives all three.  On Return Variables not
  40.      requested by JOB are not used.
  41.  
  42.      AAAA contains the upper triangle of the inverse of
  43.         the original matrix.  The strict lower triangle
  44.         is never referenced.
  45.  
  46.      DDDDEEEETTTT REAL(2)
  47.         determinant of original matrix.
  48.         Determinant = DET(1) * 10.0**DET(2)
  49.         with 1.0 .LE. ABS(DET(1)) .LT. 10.0
  50.         or DET(1) = 0.0.
  51.  
  52.      IIIINNNNEEEERRRRTTTT INTEGER(3)
  53.         the inertia of the original matrix.
  54.         INERT(1)  =  number of positive eigenvalues.
  55.         INERT(2)  =  number of negative eigenvalues.
  56.         INERT(3)  =  number of zero eigenvalues.  Error Condition
  57.  
  58.      AAAA division by zero may occur if the inverse is requested and  SSICO  has
  59.      set RCOND .EQ. 0.0 or  SSIFA  has set  INFO .NE. 0 .  LINPACK.  This
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSSSSSIIIIDDDDIIII((((3333FFFF))))                                                            SSSSSSSSIIIIDDDDIIII((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      version dated 08/14/78 .  James Bunch, Univ. Calif. San Diego, Argonne
  75.      Nat. Lab Subroutines and Functions BLAS SAXPY,SCOPY,SDOT,SSWAP Fortran
  76.      ABS,IABS,MOD
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.